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ABSTRACT 

We discuss physical experiments achievable via the monitoring of stellar dynamics near the massive black 
hole at the Galactic center with a diffraction-limited, next generation, extremely large telescope (ELT). Given 
the likely observational capabilities of an ELT and what is currently known about the stellar environment at the 
Galactic Center, we synthesize plausible samples of stellar orbits around the black hole. We use the Markov 
Chain Monte Carlo method to evaluate the constraints that the monitoring of these orbits will place on the 
matter content within the dynamical sphere of influence of the black hole. We express our results as functions 
of the number N of stars with detectable orbital motions and the astrometric precision 69 and spectroscopic 
precision 6v at which the stellar proper motions and radial velocities are monitored. Our results are easily scaled 
to different telescope sizes and precisions. For N = 100, 66 = 0.5 mas, and 6v = 10 km s -1 — a conservative 
estimate of the capabilities of a 30 meter telescope — we find that if the extended matter distribution enclosed 
by the orbits at 0.01 pc has a mass greater than ~ 10 3 M Q , it will produce measurable deviations from Keplerian 
motion. Thus, if the concentration of dark matter at the Galactic Center matches theoretical predictions, its 
influence on the orbits will be detectable. We also estimate the constraints that will be placed on the mass of 
the black hole and on the distance to the Galactic Center, and find that both will be measured to better than 
~ 0.1%. We discuss the significance of knowing the distance to within a few parsecs and the importance of 
this parameter for understanding the structure of the Galaxy. We demonstrate that the lowest-order relativistic 
effects, such as the prograde precession, will be detectable if 69 < 0.5 mas. Barring the favorable discovery 
of a star on a highly compact, eccentric orbit, the higher-order effects, including the frame dragging due to 
the spin of the black hole, will require 69 < 0.05 mas. Finally, we calculate the rate at which monitored stars 
experience detectable nearby encounters with background stars. The encounters probe the mass function of 
stellar remnants that accumulate near the black hole. We find that ~ 30 such encounters will be detected over 
a ten year baseline for 69 = 0.5 mas. 

Subject headings: astrometry — black hole physics — Galaxy: center — Galaxy: kinematics and dynamics — 
infrared: stars 



1. INTRODUCTION 

Observational programs with ten meter class telescopes, 
including the W. M. Keck Observatory and the Very Large 
Telescope (VLT), have yielded a wealth of information on 
the stellar content inside the sphere of infl uence of the mas- 
sive black hole at the Galactic cent er (GC : lGhez et alJ|1998 [ 
i Gezari et al.1 JZOOl iHornstein et all 12001 IFiger et aljl2003t 



3l et al Jl2003anGheT et al Jl2003bt ISchodel et al J l2i 
The black hole is located at the center of a compact stel- 
lar cluster that has b een the target of observational sur- 
veys for a decade (e.g.,|Krabbe et ali i 1 9951 IFiger et al] l2000t 
IGezari et al.l 1200 2). Near-infrared monitoring with speckle 
and adaptive optics techniques has recently enabled com- 
plete orbital re construction of several stellar sources o rbiting 
the black hole (lEckart et al.ll200l ISchodel et alEooH 120031 
iGhez etalJl20 03b). Sources have been monitored with astro- 
metric errors of a few milli-arcseconds ( Ghe z et al . 2003a; 
Scho del et al]l2003l). and radial velocity e rrors < 50 km s" 1 
( Eisenhau er et alJl2003l [Ghez et al.l l2003al) . allowing the de- 
tection of the accelerated proper motions of ~ 10 stars. One 
of these stars has an o rbital period of only ~ 15 yr iGhez et al.l 
2003bt ISchodel et all2003l) . 

The presence of a dark mass at the center of the Galaxy 
could in principle be inferred from the static nature of the ra- 
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dio source Sgr A* located at the center o f the stellar cluster 
jBacker & SrameUl99lllReid et alll999l) . Nevertheless, it is 
the stars with the shortest orbital periods that have provided 
unequivocal proof of the existence of a massive black hole 
and a measurement of its m ass of ~ 4 x 1O 6 M JGhez et alJ 
2003b; Schodel et al. 2003). Since, for a fixed angular scale, 

3/2 —1/2 

the orbital periods are proportional to /?q M bh ' and the ra- 
dial velocities are proportional to R^^M^ 2 where Rq is the 
heliocentric distance to the black hole and Mbh is its mass, 
the two paramet ers are no t degenerate and can be determined 
independently (Eisenhau er et alJ2003l) . 

In spite of the quality of elementary data available about 
the black hole and the bright stellar sources, the matter 
content in the vicinity of the black hole remains unknown. 
The observed stellar sources may represent only a fraction 
of the total matter content. Since the radial diffusion time 
~ 10 8-9 yr is shorter than the age of the bulge, a large number 
of massive compact remnants (5- 10M Q black holes) could 
have segregated into, and may dominate the matter density 
inside the dynamical sphere of influence of the black hole 
( Morri Jl993l: iMiralda-E scude & GoulJ2000l) . Furthermore, 
adiabatic growth of the massive black hole could have com- 
pressed a pre-ex is ting di str ibution of cold dark matter ( CDM) 
(Ipser & Sikivie Il987t lOuinlan. Hernquist. & Sigurdsson 
1 19951 jGondolo & Sffld Il999l) and stars JPeeblesI Il972t 
Young ll980l) into a dense "spike" . A variety of dynamical 
processes, however, are capabl e of de s troying such a spike 
(lUllio. Zhao. & Kamionkowskil 1200 ll iMerritt et alJ l2002t 



2 



iGnedin & Primackll2003l lMerrittll2003h . A sustained CDM 
spike would have implications for the detection of annihi- 
lation radiation for the CDM models in which the CDM 
consists of weakly-interacting massive particles (WIMPs). 
The most complete catalo gue of stars in th e central par- 



sees was compiled by Genzel et al. (2000) and Schodel et al 
20031) . In a survey of the stellar sources. iGenzel etal 



infer a spatial number density of n(r) oc r~ lA over 
the radial range 0.004 < r < 0.4 pc. Their sample was 50% 
complete for stars brighter than K ~ 18, where completeness 
is defined as the percentage of stars in the field of view that 
are detectable and thus included in the sample. Expressed in 
terms of stellar mass, the sample is 50% complete for masses 
m > 3Mq, assuming s tars on the main sequence, a distance 
to the G C of 8.0 kpc dReidlll993l) an d K-band extinction of 
3.3 mag ( iRieke. Rie ke. & PauD ll989l) . A picture is emerg- 
ing in which the brightest stars in the Central Cluster (< 0.03 
pc) are young, main-sequence stars with apparent magnitudes 
K > 13 and masses 10- 15M Q . The stars outside 0.03 pc ap- 
pear to be spectroscopically and kinematically distinct. They 
span a larger range of magnitudes K > 10 and contain ~ 40 
mass-losing Wolf-Rayet stars (e.g.. IGenzel et alJl2003al and 
R. Genzel, private communication). Unlike the Central Clus- 
ter, th ese stars appear to belong to twin, m isaligned stellar 
disks (Levin & Beloborodovl2003t IGenzel et alJ2003al) . 

The formation of the observed young stars with x 100 
larger specific binding energies relative to the black hole 
than that of th e nearest observed a ccumulation of molec- 
ular gas (e.g., LTackson et al.l 119931) presents a challenge 
to st ar formation theories and is a persistent puzzl e 
(e.g., lMorrislll993l iGhez et alJl2003at IGenzel et al.ll2003al) . 
A number of mechanisms for the formation and migra- 
tion of stars in the tidal field of the massive black hole 
have been proposed (lGerhardl l2001l iGould & OuilleijpQQ^I 
Hansen & Milosavlievic 2003; Levin & Beloborodov 2003; 
iKim & Morrisl2003HMilosavlievic & Loebl2004l) . mile the 
mechanisms have important implications, they are also each 
deficient in at least one way. 

The re is a dear th of giants in the GC region ( Eck art etail 
1995). Recently, iFi ger etalJ (120031) measured the radial ve- 
locities of 85 cool, normal giant stars with projected dis- 
tances from the central region between 0.1-1 pc. They find 
nearly complete deficiency of giants with large radial veloci- 
ties (V rac i > 200 km s -1 ). Since a star in a circular orbit at a dis- 
tance of 0. 1 pc from the black hole has velocity ~ 400 km s~ , 
the absence of any such stars with comparable radial veloci- 
ties indicates that the observed giants are indeed limited to the 
region outside the central ~ 0.5 pc. 

While the measured stellar density profile of the Galac- 
tic bul ge is consisten t with that of a singular isothermal 
sphere (Becklinl 119681) . the profile in the central parsec is 
not well known, especially for the lower-mass stellar popu- 
lations. Assumiiigj£k«atiOTthat is driven by two-body pro- 
cesses, Bahcall & Wolf ( 1976) showed that the equilibrium 
phase space distribution for a population of equal mass stars 
is a power law in density p oc r" 7 / 4 . For a multimass distri- 
bution the lighter stars are less centrally concentrated, result- 
ing in a power-law profile that ranges from r~ 3 ' 2 for the least 
massive species to r~ 7 / 4 for the most (|Bahcall & Wolfl l 19771 
Murp h*v~Cohn. & Durisenl 19911) . A coeval family of stars in 
the central region has reached equilibrium only if it is older 
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where a is the local linear stellar velocity dispersion, /«* is the 
mass of a typical field star, p is the local stellar density, and 
In A is the Coulomb logarithm. 

Since the main sequence lifetime of stars more massive 
than ~ 2Mq is shorter than fg, young massive stars in the 
GC are not relaxed; their distribution is primarily a reflec- 
tion of their formative conditions. While lower mass dwarf 
stars are sufficiently old to be relaxed in the central potential, 
their distribution in the innermost region could be affected 
by an abundance of stellar mass black holes (5 - IOMq). As 
products of normal stellar evolution, stellar mass black holes 
sink in the potential of the ma ssive black hole (Morrisl [l993l: 
Miralda-Escude & Gould 2000) and displace the less massive 
stars and remnants. 

Speckle imaging and more recently adaptive optics with the 
Keck and VLT have provided several milliarcsecond astrome- 
try, enabling the detection of proper motions within the inner 
0.5 pc and accelerated proper motions of ~ 10 stars within 
the inner 0.05 pc. Radial velocities with spectroscopic preci- 
sions of Sv ~ 30 km s -1 have also been obtained for the star 
S0-2, which has been monitored for over 70% of its orbit in- 
cluding pericenter passage at ~ 130 AU from the black hole. 
These observations have enabled the black hole mass a nd GC 
dista nce to be meas ured to within ~ 10% (Ghez et al. 2003b; 
Schodel et alJ200l . 

Here we examine the extent to which one can probe the 
GC potential by monitoring stars with a diffraction-limited, 
next generation, extremely large telescope (ELT). As com- 
pared with current 10 m class telescopes, the finer angular res- 
olution of an ELT enables the orbital motions of many more 
stars to be detected, each at greater astrometric precision, 89, 
and spectroscopic precision, i5v. Given the range of possible 
sizes of future telescope and given the uncertainties in the ul- 
timate capabilities of a specific telescope class (e.g., 30 meter 
telescopes) we choose to express our results not as functions 
of the ELT aperture but rather as functions of 89, Sv, and the 
number TV of stars with detectable orbital motions. We take 
69 = 0.5 mas and 6v = 10 km s -1 as our fiducial model, cor- 
responding to a conservative estimate of the capabilities of a 
telescope with a D = 30 m aperture. We show that TV scales 
with telescope aperture as TV ~ 100(Z)/30 m) 2 . We demon- 
strate that with an ELT one can measure the density profile 
of a dark matter spike and those general relativistic effects 
that scale as (v/c) 2 , where v is the speed of a star and c is 
the speed of light. Furthermore, we show that the distance to 
the GC will be measured to remarkable precision. This will 
help place tight constraints on models of the overall Galac- 
tic structure. We also show that with an ELT one can de- 
tect the gravitational interactions between monitored stars and 
the background massive stellar remnants that accumulate near 
the central black hole. Such interactions may probe the mass 
function of the stellar mass black holes thought to dominate 
the matter density in the region. 

The paper is organized as follows. In §[2] we calculate the 
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number of stars with accelerated proper motions that can be 
monitored with a given ELT based on its astrometric, spec- 
troscopic, and confusion limits. We also describe a realistic 
monitoring program and demonstrate that confusion with the 
infrared emission from Sgr A* is unlikely to affect an ELT's 
ability to measure stellar motions. In §[3]we model the orbital 
data and estimate the magnitude of various non-Keplerian ef- 
fects including Newtonian retrograde precession due to ex- 
tended matter, relativistic prograde precession, precession in- 
duced by the coupling of orbits to the spin of the black hole, 
and the Roemer time delay. In § 13.61 we consider the effect 
of stellar interactions on the motion of the monitored stars. 
Specifically, we estimate the rate at which discrete stellar en- 
counters result in detectable changes of orbital motions. In §|4] 
we discuss a method for generating mock ELT orbital data and 
describe a computational technique for estimating uncertain- 
ties in the orbital parameters. The results of our calculations 
are given in §|5] Finally, in §[6] we discuss astrophysical ap- 
plications of the proposed observations. 

2. OBSERVING STARS IN THE CENTRAL ARCSECOND WITH AN 

ELT 

The purpose of this section is to estimate the number and 
distribution of stars whose orbital motions can be detected 
with an ELT and to determine the astrometric and spectro- 
scopic precision to which their motions can be measured 
(§ 12. U . The latter are determined by the specifications of 
the telescope and the properties of the stellar population at 
the GC. Several factors complicate the monitoring of orbits 
within the central arcsecond. The greatest obstacle to detect- 
ing and following hitherto unseen stars is stellar crowding. 
Light contamination from nearby bright stars as well as the 
light from underlying faint stars flood the pixel e leme nts and 
impose a limit to the faintest detectable star. In § 12. 21 we esti- 
mate the minimum luminosity permitted by the crowding and 
thereby obtain an estimate of the number count of stars with 
observable orbital motions. 

2.1. Astrometric and Spectroscopic Limit 

With adaptive optics, an ELT will operate near its diffrac- 
tion limit in the /f-band. By determining the centroid of im- 
ages, the measured relative positions of stars are a factor of 
~ 20-40 more definite than the images' diffraction limit. For 
instance, the diffraction limit of Keck is ~ 50 mas while the 
astrometric error of a bright star near the GC as seen by Keck 
is ~ 1-2 mas. Naively, the expected astrometric limit of an 
ELT with D = 30 m is therefore Sdjo ~ 0.5 mas. 

In practice, the astrometric limit achievable with adaptive 
optics depends on whether atmospheric fluctuations or cen- 
troid measurement errors dominate the signal. At the GC 
the separation between the guide star needed for the adap- 
tive optics infrared wavefront sensor and the star under study 
is typically ~ 5", corresponding to a separation of 0.25 m at 
the top (~ 10 km) of the atmosphere. As long as the tele- 
scope aperture is larger than this separation, as is the case for 
Keck and an ELT, the atmosphere dominates and the astro- 
metric precision scales with the telescope diameter as Z) 2 / 3 
iShao & Colavitall 992). A 30 meter ELT is expected to have 
an astrometric limit that is 3 2//3 w 2 times smaller than Keek's 
for K < 24. We therefore adopt 89^ ~ 0.5 mas in our cal- 
culations, though we consider this a conservative estimate; a 
30 meter ELT may attain an astrometric limit as small as 0. 1 



mas. 3 

With an adapt ive-optics-fed spectrometer on Keck, 

iGhez et all i2003al) detected spectral absorption lines in 
the star SO-2 at a spectral resolution ofi? = A/AA~ 4000, 
yielding a radial velocity measurement with an error of 
40kms _1 (see also Eisenhau er et alJl2003t) . Integral field 
spectroscopy in the near-IR with a 30 meter ELT is expected 
to enable measurements with R ~ 1 -2 x 10 4 , suggesting that 
velocity errors of 5v3o ~ 10 km s -1 are attainable. This too 
is a conservative estimate as a 30 meter ELT may achieve 
velocity errors more than an order of magnitude smaller (D. 
Figer, private communication). As we discuss below, an 
ELT will be able to detect stars that are fainter and hence 
cooler than those currently detectable. Cool stars exhibit 
rich spectral features including possible molecular lines, 
enabli ng high s pectral resolution studies. For example, 
Figer et al. (2003) obtained radial velocities for 85 cool stars 
in the central parsec of the Galaxy with velocity errors of 
- 1 kms" 1 . 

Although an ELT's astrometric and spectroscop ic lim its 
may differ from the above estimates, we show in § 14.21 that 
the uncertainties in the model parameters extracted from the 
monitoring data, such as the distance to the GC and the ex- 
tended matter profile, scale almost linearly with the measure- 
ment errors. The constraints on the parameters for different 
values of 86 and 8v can therefore be readily inferred from our 
results. 



2.2. Confusion Limit 

The brighter stars wash out the signal of fainter stars, 
thereby limiting the luminosity of the faintest observable star. 
This limit depends on the telescope optics (e.g., angular res- 
olution) and on the stellar luminosity function (LF). Using 
measurements of stellar photometry near the GC, we now es- 
timate the minimum luminosity that a star at the GC can have 
and still be identified and monitored with an ELT. For a given 
star of luminosity / and for a given /if -band stellar LF, we de- 
termine the integrated flux from all nearby background stars 
with luminosity < /. At some minimum luminosity, the emis- 
sion from a single star is comparable to the background emis- 
sion; this lu minosity sets the confu sion limit. 

Following Takeuchi et al] d2001h and references therein, let 
Xs = Sh(0,4>) be the response of the telescope to a source of 
flux density S at an angular position (8, </>) from the line-of- 
sight axis to the center of the source. h(9,<fr) is the point- 
spread function (PSF) of the telescope, normalized to unity at 
the center. Since all sources at the GC are essentially at the 
same distance, we can instead express the response in terms 
of stellar luminosity I, i.e., let x = lh(9,tfi). The variance in 
the telescope response due to crowding is the confusion noise 
a. To detect a source with high statistical significance, its 
luminosity must be greater than some cutoff l c , or equiva- 
lently, x must be greater than a response cutoff x c . Defining 
q = x c /a, a source is above the confusion limit if its signal-to- 
noise S = x/ a > q, where we take q = 5. 

If the number of stars per square arcsec with luminosity in 
the range (l,l+dl) is dN= a§(l)dl, where a is the normaliza- 
tion of the LF and $(/) is its shape, then the mean number of 

3 For example, see http://tmt/ucolick.org/reports_and_notes/index/htm 
Report No. 34. 
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FIG. 1 . — The 7Y-band magnitude limit and number of stars with detectable orbital m otions as a fun ction of the aperture of a diffraction limited ELT. Results 
are shown for powe r-law 7?-band luminosity functions normalized to observations by Schodel et al. 1 2003 ) with slopes matching the ~ 2cr range found by 
IGenzel et all f2003a1. The uncertainty in the number counts is considerably smaller than the uncertainty in the magnitude limits. 
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where the integral is over the solid angle of the PSF. The con- 
fusion noise a due to all sources fainter than x c is then 



a 2 = 



x R(x)dx. 



(3) 



Since we are interested in calculating the cutoff response 
of a given detector for a given LF, we need to solve for the 
confusion noise. Assuming a power law LF of the form 
dN/dl = = al~ n we have 
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In this paper we only consider power-law LFs, though one 
can obtain an expression for a for general forms of the LF 
JFranceschini et alH"9 89). 

For a Gaussian PSF h{9,f) = h(8) = exp[-(4ln2)((9/6>o) 2 ], 
where 6>o is the PSF's full-width at half-maximum. This gives 
^eff = ?r(9g/ [(41n2)(rj- 1)] and thus 



3 _„ x i/(n-D 



irdia 



(41n2)(77- 1) 



1/(17-1) 



(7) 



We now estimate the value of a for the Keck and an ELT. 
In the /T-band, 9o — 50 mas for Keck and 6q — 15 mas for a 
D = 30 m ELT. We also ne ed the A"-band lumino sity function 
(KLF) of stars at the GC. IGenzel et all J2003al) find that the 
KLF within 1.5" of the GC is well described by a power- 
law with slope (3 = dlogN/dK = 0.21 ± 0.03 where K is the 
apparent magnitude in the /T-band. We consider KLFs with 
slopes within the ~ 2<r range 0.15 < j3 < 0.27, which in terms 
o f 77 = 1 + /3/0.4 c orresp onds to the range 1.38 < i] < 1.68. 

ISchodel et alJ J2003I) measured the photometry of more 
than 40 stars in the central arcsec, 29 of which reside within 
0.8" (~ 6000 AU). We normalize the KLF to these 29 stars. 
We limit our analysis to these innermost stars since the KLF 
inside 0.8" appears to differ from that outside this region (see 
§Q. We do not attempt to account for a possible radial de- 
pendence but instead assume the KLF is constant. 

Of the 29 stars, the brightest has apparent magnitude K = 
13.4 and the di mmest K= 17.3. Assuming a /T-band extinc- 
tion of 3.3 mag jRieke. Rieke. & Paul 1 989) and a distance to 
the GC Rq = 8 kpc, these apparent magnitudes correspond to 
T^-band luminosities of / m ; n = O.8L0 and l max = 28L©. For a 

given 77 we calculate a = N obs (l - rj) / '(Z^2 - ) where A^ obs = 
29/7r(0.8") 2 and by equation @ solve for a. Integrating the 
luminosity function over stars brighter than x c = qa yields the 
number count of detectable stars N(l > l c ) ~ ax 1 ^ /{rj- 1). 

In Figure ^ we show how the /T-band magnitude limit 
and number of stars with detectable orbital motions (those 
within 3000 AU of the GC; see below) scale with the aperture 
of a diffraction limited ELT assuming (3 = 0.21 ± 0.06 (77 = 
1 .53 ± 0.15). Since by equation /V oc xl~ n oc 8^ 2 , we find 
that N ~ 100(D/30 m) 2 . Furthermore, because x c oc a 1/(ri ~ l) , 
N(l > l c ) is not very sensitive to the value of x c for a fixed 77. 
Therefore, when the above analysis is performed on a subset 
of the 29 stars within 0.8" (e.g., stars within 0.4" or alterna- 
tively stars brighter than K = 16), the derived number counts, 
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unlike the magnitude limits, do not change significantly. The 
number counts we derive for an ELT are therefore robust even 
though the magnitude limits are subject to some uncertainty. 

To extract orbital parameters the acceleration of a star in the 
plane of the sky must be detected, i.e., it must be greater than 
the threshold acceleration For Keck ~ 1 -2 mas yr" 2 
while for a 30 meter ELT £, ~ 0.5 mas yr~ 2 . The acceler- 
ated proper motion is detectable over the entire orbit if the 
acceleration at apocenter exceeds the threshold. For a face 
on orbit this requires a(l +e) < (GM/&R0) 1 ' 2 . Thus the ac- 
celeration will be detectable with a 30 meter ELT over the 
entire orbit if a < 3000 AU (period < 80 yr). To construct our 
mock stellar orbits to simulate observations that can be made 
with an ELT, we only consider orbits satisfying this constraint. 
As Figure[2shows, within 3000 AU, approximately 100 stars 
are brighter than a 30 meter ELT's confusion limit. Further- 
more, since the surface density of stars is ~ 200 arcsec -2 , one 
does not expect to observe stars with apocenters smaller than 
- 300 AU (~ 0.04") with such an ELT. We therefore con- 
clude that a 30 meter ELT will detect the accelerated motion 
of ~ 100 stars with semi-major axes < 3000 AU and apocen- 
ter distances > 300 AU. 

Another related issue is the frequency with which an ELT 
will measure positions and radial velocities for the N mon- 
itored stars, given a reasonable commissioning of ~ 10 GC 
exposures per year. As it will be equipped with an integral 
field unit spectrometer an ELT can obtain simultaneous spec- 
tral and spatial data over a relatively large region of sky. It 
is possible for it to measure positions and velocities for all N 
stars in a single image. This suggests that a dedicated observ- 
ing program can reasonably obtain ten measurements per star 
per year. 

At the high levels of precision obtainable with an ELT, 
orbital parameter constraints should scale with the measure- 
ment errors a (i.e., 89 and 6v) and the number of stars N as 
a/N 1 ' 2 . Based on the above discussion we therefore expect 
the parameter constraints to scale with telescope aperture as 
~ D~ 2 / 3 /D ~ ZT 5 / 3 . We verify this relation in our numerical 
simulation results described in §|5] 

2.3. Central Point Source — Sgr A* 

At rad io wavelengths Sgr A* is detect e d as a non- 
therm al dBeckert et alJ 119961 ISerabvn et all 119971). com- 
pact (iRogers et alJ 119941) , sta tic JBacker & Sramekl |1999t 
iReidetalJ 119991) . variable JZhao. Bower. & Gossl 120011) 
source. An X- ray source coinciden t with Sgr A* has also 
been detected (Bag anoff et al.l 120031) and consists of a re- 
solved, steady-state component with size ~ 1" and an unre- 
solved flaring component that increases in flux density by an 
order of magnit ude over the course of a few hours roughly 
once per day ([B aganoff et al.l 1200 It |Qoldwurm et alJ 120031 
Porquet et all2003l) . 

Recently, a near-infrared counterpart to Sgr A* , located 
within a few mas of the dynamically determined black 
hole position, has b een observed in the H ba nd (1.7fzm; 
[ Genzel et all 12003b!) and the U band (3.8/im; iGhez et all 
2004). Like the X-ray emission, the infrared emission con- 
sists of a quiescent component and a variable component. The 
latter exhibits flux densities that increase by a factor of a few 
over the course of tens of minutes to one w eek, with possi- 
ble signs of periodicity (Gen zel et al.ll2003 li). The observed 
L'-magnitudes are in the range 12.2- 13.8, corresponding to 
L'-luminosities ~ 10- lOOL©. 

Although the current sample of stars with detected accel- 



erated motion are brighter than the Sgr A* infrared emission, 
the stars detectable with an ELT will have comparable lumi- 
nosities. A star that passes near the black hole can therefore 
be confused with the emission of Sgr A* . Conservatively, 
such confusion limits monitoring when the projected separa- 
tion between a star and Sgr A* is smaller than the resolution 
of the detector. For a 30 meter ELT operating at the diffraction 
limit in the K band, this co rresp onds to ~ 15 mas (120 AU). 
However, as we found in § 12.21 confusion with nearby stars 
precludes such a telescope from detecting orbits with apocen- 
ters smaller than 300 AU (~ 40 mas) and most of the moni- 
tored stars do not therefore pass within 15 mas of Sgr A* . Of 
those that do, most spend only a small fraction of their total 
orbital period that close to the black hole; e.g., a star with a 
semi-major axis of 200 AU and eccentricity 0.9 is within 15 
mas of the black hole for only 10% of its orbital period. The 
same arguments hold for other ELT apertures. Therefore, the 
infrared emission from the black hole will not significantly 
impair orbital monitoring with an ELT. 

3. ORBITAL DYNAMICS 

While current observations of stellar proper motions near 
the black hole at the GC are consistent with motion around 
a Newtonian point mass, we show that with an ELT non- 
Keplerian motions are going to be detectable. There are var- 
ious effects that cause deviations from Keplerian motion, in- 
cluding the Newtonian retrograde precession (NRP) of an or- 
bit due to the presence of an extended matter distrib ution 
(§ 13. 21 . the relativistic prograde precession (RPP; § 13. St - 
an d th e frame dragging effects related to the black hole spin 
(§ !3.4t . In addition, we account for an apparent d efor mation 
of the observed proper motion ("Roemer effect"; § 13.51 ) due to 
the differences in light travel times at different locations along 
the orbit. A discussion of the effects of encounters betw een 
monitored stars and background stars is given in § 13.61 We 
now describe the orbital equations of motion and estimate the 
magnitude of the various non-Keplerian effects. A number of 
relativistic effects, including those we consider below, are dis- 
cussed in Pfahl & Loeb (2003) in connection with long-term 
timing observations of a radio pulsar that might be detected in 
a < 100 year orbit about the GC. 

3.1. Equations of Motion 

We found in §[2] that we do not expect a 30 meter ELT to 
detect orbits with apocenter smaller than ~ 300 AU due to 
confusion noise. Assuming orbits uniformly distributed in <? 2 , 
the probability that a given star has e > 0.99 is 2%. Since 
most of the ^100 stars such an ELT monitors will have semi- 
major axes > 1000 AU, it is unlikely that any will have peri- 
center distance smaller than a few AU. As a result, the ratio 
of the Schwarzschild radius to the pericenter distance of the 
stars will satisfy R s /r p < 0.05, or expressed in terms of the 
stellar velocity at pericenter, v p /c< 0.2. The post-Newtonian 
approximation to the geodesic equations that is accurate to or- 
der (v/c) 2 provides an adequate description of the stellar or- 
bits given the observational precision expected with an ELT. 

The geodesic equation for test particles orbiting a spheri- 
cally syjnm^tricji^sj^^ approxima- 
tion jWeinbergll972HRubilar & Eckari2001l) . 



- [2V$ 2 + v 2 V$-4v(v ■ V)$-v x (V x 0] /c 2 , (8) 
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where v = dx/dt is the velocity vector, $ is a time- 
independent gravitational potential, and ( = 2G(x x J) /r 3 is 
a vector potential associated with the spin J of the gravitat- 
ing mass, which we assume is constant with time. We as- 
sume the density distribution of the extended matter at radii 
traversed by the stars and smaller is a power-law profile 
p (r) = poir/rof 1 . Input model parameters are described in 

The relativistic effects include corrections to the orbital dy- 
namics and to the observed motion due to propagation effects. 
The former class includes the RPP and frame dragging while 
the later class includes the lowest order (v/c) Roemer time 
delay and such (v/c) 2 effects as time dilation, gravitational 
redshift, and the Shapiro time delay. Since the (v/c) 2 propaga- 
tion effects each have different functional dependences on the 
orbital parameters (see e.g., iPfahl & Loebl l2003). including 
them may break degeneracies, though they may also weaken 
the sensitivity to some parameters. However, our interest in 
the (v/c) 2 relativistic effects is primarily connected with the 
ability of an ELT to probe general relativity on the scale of 
a massive black hole rather than with the effects' potential 
utility for parameter estimation. Since v p /c < 0.2, relativis- 
tic effects help to constrain the orbital parameters by at most 
a few percent. We therefore chose not to include the (v/c) 2 
propagation effects in our analyses. The analysis of actual 
data obtained from an ELT must, however, account for all the 
relativistic effects. 

3.2. Newtonian Retrograde Precession 

The NRP was discus sed in the context of the GC by 
Rubilar & Eckart (2001). An extended matter distribution 
causes stellar orbits to precess due to differences in the 
amount of mass that is contained between the apocenter and 
the pericenter radii. In the numerical calculations that fol- 
low, we determine how much an orbit precesses due to the 
extended matter by solving equation (|8j. Here, however, we 
obtain an estimate of the magnitude of the precession for stars 
at the GC by considering the potential of the extended matter 
to be a small correction (5$ to the potential of the black hole. 
Expanding the tota l potential to linear order in i5$, the angular 
shift per period is ( Landa u~& Lif shitzl 19 60) 

A ^=]lz(zl r 2 5<S>d^j, (9) 

where L is the orbital angular momentum per unit mass, 
r = a(l -e 2 )/(l +ecostp), and ip is the phase of the orbit. 

If 7 < 2 we have (5$ = /3r 2 ^, where ft = GM ext /(2- 7 )ro~ 7 
is a constant and M ext = M ext (r < ro) is the extended mat- 
ter mass within rp. The o rbital precession is then given by 
JMunvaneza & Viollierll99"9l) 

A0 New t= kl-e 2 )]%( 7 ,e), (10) 



GM bh 



where 



and 



g( 7 ,e)= —^-(4-7) [/ 4 - 7 (e)-/ 5 - 7 (e)] 
+(7-2 7 )/ 4 - 7 (c) 



In(e) : 



d^p 



(11) 



(12) 



/ () ( 1 + e cos ip) n 

Assume that the extended matter consists of stars with 
7 = 7/4 and pi pc = 2 x 10 5 M Q pc -3 (see § Consider an 



SO-2 like orbit with a semi-major axis of 0.005 pc and eccen- 
tricity e = 0.9. The enclosed stellar mass at apocenter and 
pericenter are 6000M Q and 15OM . Solving equation dlOt 
yields a precession per revolution of A^Newt ~ 0.08°, corre- 
sponding to an apparent angular apocenter shift of roughly 
A<p^ ewt a(l +e)/Ro « 0.3 mas. Thus, a few SO-2 like orbits 
with astrometric errors of 0.5 mas provide a meaningful con- 
straint on the stellar distribution within the inner few milli- 
parsecs. If the density of the dark matter cusp at the stellar 
positions exceeds ~ 1O 8 M pc" 3 , then it too will produce a 
detectable precession; it will not be ea sily distinguished from 
the stellar contribution (however see § 13.61 . 



3.3. Relativistic Prograde Precession 

T he RPP causes a pericenter advance per revolution of 
(see IWeinberglll972[) A0 plo = 3nR s /a(l-e 2 ), where R s = 
2GMbh/c 2 is the Schwarzschild radius of the black hole. The 
magnitude of the effect is oc (v/c) 2 . The apparent apocenter 
shift per revolution caused by the RPP is As A<j> pm a(l + 
e)/Ro = 3ttR s /Rq(1 - e), which corresponds to an apparent 
shift of ~ 1 mas for the star SO-2. Although the RPP h as an 
additional factor of (v/c) relative to the Roemer effect (§ 13. 51 . 
this attenuation can be compensated by having a few high ec- 
centricity stars in the sample. Furthermore, unlike the Roemer 
effect, the RPP shift is to first order independent of the semi- 
major axis and is therefore equally sensitive to stars at all radii 
(although stars at large radii also have long periods). Consider 
an orbit seen face on and observed for N or \, complete periods. 
Since the precession angles per revolution add linearly, the 
signal-to-noise from the RPP is S pro ~ AsN OI b /S9, or 



J pio 



^0.1 



1-e V4x 10 6 M P 



Ro 



8kpc 



1 mas 



(13) 



In a sample of 100 stars observed with astrometric errors of 
0.5 mas and having an eccentricity distribution uniform over 
e 2 , we expect on average eight stars with e > 0.96. If only one 
such a star is followed over just a single period, the RPP shift 
will be measured to 5-a accuracy. 

3.4. Frame Dragging 

For a spinning black hole, frame dragging effects also cause 
a precession of the pericenter. The spin precession per revolu- 
tion for a star orbiting a black hole with spin angular momen- 
tum J is given approximately by (see Weinberg 1972, equa- 
tion (9.5.22); note different notation) 



spin ' 



-8tt; 



GM bh \ 



cL 

2V2jA0 pro 



cos?/; 



R s 
a(l-e 2 ) 



cos-0, 



(14) 



where ijj is the angle between the orbital angular momentum 
vector and the black hole spin axis and < j = cJ /GM 2 ^ < 1 
is the black hole spin parameter. 

The black hole spin induces an apocenter shift that is 
smaller than the RPP shift by a factor of ~ v/c. Even if 
the black hole is maximally spinning (j =1), the shift rep- 
resents only a 5% contribution on top of the RPP for a star 
with a = 200 AU and eccentricity e = 0.92. For an orbit ob- 
served face-on the signal-to-noise from a spin-induced apoc- 
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enter shift is 



2ir\/2j 



R s \ V2 N orb 



COS-0 



i 0.001 



V(l+e)(l-e) 3 V4 x 10 6 A/, 



3/2 



( " ) 

V1000 AU/ 



1/2 / /? 



8kpc 



56* 
1 mas 



■ (15) 



For example, a 5-er detection is achieved with an ELT with 
59 = 0.5 mas if a star with a = 300 AU and e = 0.99 is moni- 
tored for three complete orbits. We expect a 30 meter ELT to 
detect one star with a semi-major axis that small (§ As- 
suming eccentricities uniformly distributed in e 2 the proba- 
bility that star has e > 0.99 is only — 2%. If 59 = 0.05 mas, a 
star with a = 300 AU and e = 0.95 will yield a 5-er detection 
after being monitored for three complete orbits. Since so high 
a resolution requires an ELT with aperture D ~ 100 m, there 
will be several stars with a < 300 AU and of these ~ 10% will 
have e > 0.95. Detecting such a stellar orbit is therefore not 
unlikely. The spin-induced orbital precession thus requires an 
astrometric precession of 59 < 05 mas (see also Ja roszvnskH 
1998: lFragile & MathewsKOOTO . 

3.5. The Roemer Time Delay 

For orbits with non-zero inclination, the distance between 
the Earth and star, and hence the difference in time between 
stellar emission and observation, varies with orbital phase. 
This time delay, given by Af = f D b s - t em = z(t em ) / c, where z(t) 
is the relative distance between the star and the massive black 
hole, was first recognized by Roemer in 1676 in application to 
the phases of Jupiter's moons. Unlike the relativistic Doppler 
effect which includes corrections of order (v/c) 2 and higher, 
the Roemer delay is the classical Dopp ler effect which only 
includes terms up to order v/c (see, e.g., Loeb 2003). The de- 
lay has a magnitude corresponding to a few percent of a year 
for an SO-2 like orbit, and is observed as an additional shift 
in the apparent stellar position with time, As(f). For a cir- 
cular orbit seen edge-on the stellar positions z(t) and s(t) are 
sinusoidal so that 

As(t)/a = cos(wf b s )-cos(u;fem) 

= cosMfem + a/c sin(wf em ))] - cos(wf em ) 



WO . 2, \ 

sin (ut em ), 

c 



(16) 



where oj = 2tt/P and we used the fact that for orbits at the GC 
v <C c. The maximum shift, in units of the semi-major axis, is 
therefore v/c. For non-zero eccentricity and arbitrary inclina- 
tion the star's projected position and distance as a function of 
time are (see e.g., Murray & Dermott 1999) 

f*(?n r cos [if {t) + a] ) 

< y(t) \=r(t)l sin [<p(t) + a] cos i \ (17) 

I z(t) J { sin[<£>(f) + a]sin/ J 

where i is the inclination, ip(t) is the orbital-phase (i.e., the 
true anomaly), a is the argument of pericenter, and we chose 
the reference direction so that the x-axis coincides with the 
longitude of ascending node. The Roemer shift is then | Aj| = 

1 /2 

(Ax 2 + Ay 2 ) where Ax = x(t \, s ) — x(t em ) and similarly for 
Ay. To linear order in v/c the orbit-averaged Roemer shift 
can be written as 

(\As(t)\/a)~— siniV/(e,i,a) (18) 



where f(e, i, a) is a factor of order unity. For the two extreme 
cases a = and a = ir/2 (corresponding to the line-of-node 
along the major-axis and minor-axis, respectively), f(e,i,a) 
is given by 



f(a = 0) = ^-f- \le A cos 2 i (3 - VT^e 2 ' 



+3e z 



2-\/l-e 2 -cos 2 /(4-3Vl-e 2 



-6sin z / 1 



f(a = n/2)-- 



2e 4 



VT^)}, 

{-2e 4 (2-Vl^e 2 



(19) 



+e 



10-7Vl-e 2 - 



-cos 2 / (4-VT- 



-6sin z / 1 



-x/T^ 2 )}. 



(20) 



An ELT will be able to detect the effect of the Roemer delay 
in orbits at the GC. The signal-to-noise from 2V b S observations 
of an orbit measured with astrometric errors 59 is approxi- 



mately Sdelay ~ (As)N^ / R 89, 01 

Sdelay ~ 0.8vW o bs S11X I y/f(e,i,a) ( 



4x 10 6 M, 







1/2 



Rt 



1000 AU 
l 



1/2 



8kpc 



1 mas 



(21) 



If, e.g., we pick i ~ 7r/3, a = 0, and e ~ 1 /v 2, an astrometric 
error of 0.5 mas, a mean semi-major axis of 1000 AU, and 10 
observations per star, then we can detect the delay to Sdeiay ~ 5 
with roughly 10 stars. We therefore expect the Roemer delay 
to be detectable in such an ELT's sample of ~ 100 stars and 
the effect must be taken into account during parameter esti- 
mation. 

3.6. Interstellar Interactions 

In the previous sections we described the motion of a star 
in the potential of a black hole and a smooth distribution of 
extended matter, including stars, remnants, and dark matter. 
This approximation ignores the fact that the potential due to 
stars and remnants is the sum of discrete point-mass poten- 
tials and is therefore not perfectly smooth. The stars experi- 
ence perturbations due to nearby encounters with individual 
stars and due to fluctuations in the potential arising from all 
stars. These perturbations cause a star's orbital parameters 
to change with time. The magnitude and the rate of these 
changes depend on the stellar mass function since the per- 
turbations are sensitive to the characteristic mass of the field 
stars. Thus, measuring the effects of stellar encounters is a 
probe of the mass function in the central parsec. It also breaks 
the degeneracy between the contributions of stellar matter and 
dark matter to the Newtonian orbital precession. Encounters 
may also be a source of noise in measurements of orbital pa- 
rameters such as the black hole mass and distance to the GC. 
While we do not include the effects of encounters in our nu- 
merical calculations presented in §[5] we now estimate their 
magnitude and demonstrate that the encounters might be de- 
tectable with an ELT and present a powerful probe of the mass 
function of stellar remnants at the GC. 

An encounter between a test star of mass mj and a field star 
of mass m, with impact parameter b induces a change in the 
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test star velocity given by (see, e.g.. lSpitzerl 19871) 



5v = 



2m,v re i 



m, + m 



1 + 



-1/2 



(22) 



where bo = G(mi+ntj)/v 2 el and v re i is the initial relative veloc- 
ity of the stars. The encounter induces a change in the test 
star's velocity distinct from that due to orbital motion around 
the black hole. We solve for the maximum impact parameter 
£>max such that an encounter induces a change in velocity of the 
test star larger than the minimum detectable change <5v m ; n . For 
uncorrelated position measurements the minimum detectable 
change in velocity is dv m [ n ~ 86Ro/ \fN h s T , where T is the 
time baseline over which the orbit is monitored, and /V bs is 
the number of position mea surements taken in time T. As- 
suming 56 = 0.5 mas f§l2~Tl. T = 10 yr, and N obs = 100 yields 
£>Vmin ~ 0.2 km s _I . By equation (1221 we have 




2m,- v re i 



1 



2Gnii 



(23) 



a 

where the approximation assumes v re i ^> (5vmi n and mj < m\. 
Assume v re i ~ v p = (GM(1 +e)/(l —e)a) for an encounter 

1 /2 

near pericenter and v re i ~ (GM/a) for an encounter near 

apocenter. If we take (5v m ; n ~ 0.2 km s -1 and m, = 10M Q then 
for an SO-2 like orbit, b max ~ 10 AU at pericenter and b max ~ 
50 AU at apocenter. 

We ignore the effect of the black hole on the encounter and 
treat the interaction between the stars as a two-body problem. 
This is a fair approximation as long as the duration of the en- 
counter is much shorter than the time scale over which the 
orbital velocity changes significantly due to the influence of 
the black hole. At pericenter passage, where the orbital accel- 
eration is greatest, the orbital time scale is t p ~ (l-e) 3 ' 2 P, 
where P is the orbital period. The two-body approxima- 
tion is valid as long as the duration of the encounter satisfies 
?enc ~ frmax/Vrei <^ tp- F° r an SO-2 like orbit t p ~ 0.5 yr while 
by equation J23I ? enc < 0.01 yr even for m, = 20M Q . 

Next, we estimate the rate at which encount ers b < b max 
occur for a star on a given orbit (see, e.g., I Yul 120031) . Let 
Tij{r,\j,t)dmi be the rate at which a star with mass mj at po- 
sition r with velocity \j at time t encounters stars with masses 
in the range range m, — + + Assume the number den- 
sity of stars is spherically symmetric and follows a power 
law v(r) = va{r/rh)~ a . The phas e-space di stribution func- 
tion of the stars is given by (Magorri an & TremaindH999i) 
/(£) = h(a)£ a - 3 / 2 , where £ = *(r)-v 2 /2 and #(r) is the rel- 
ative gravitational potential at r, while 

r(a+l) 



A(a) = (27ro0 



2x-3/2 



V0 



IXa-1/2) 



2a+3 



(24) 



with <Ji, the linear stellar velocity dispersion outside the sphere 
of influence of the BH r/, ~ 1 pc. The rate of detectable en- 
counters in the mass bin is then 



dvi f(r,Vi) 



= 2nKh(a) 



dv re lV re iS(Vrel) 



|Vj-Vj| 



a-3/2 



- (vi + Vj-\Vi-Vj\)dVi, 



(25) 



where the cross section for detectable encounters S = 7rfr^ ax , 
andK = 4nG 2 ml/5v 2 min . 

We now determine the rate at which stars that will be mon- 
itored with an ELT undergo detectable encounters. The inte- 
gral in equation d25i is most easily evaluated in the special 
case a = 3/ 2, which is co mpatible with current observational 
constraints JGenzel et alJl2~0 03a). To obtain a rough estimate 
of the rates, consider the case a = 3/2 and assume the back- 
ground stars all have identical mass not smaller than that of 
the test star (e.g., they are a population of stellar mass black 
holes). By equation (I25> 



r/iw) = 



1 



(2*--vj), 



(26) 



and upon averaging over the orbital phase 



i r 

Tj(e,a)=- / Tj{r,yj,t)dt 
r Jo 

5-k\[2vo ( Grtii \ GMb 



2a\ 



Sv„ 



(27) 



Assume the stars monitored with an ELT have an eccentric- 
ity distribution uniform in e 2 (isotropic velocity ellipsoid) so 
that dN / de da oc ea 2 ~ a = e^fa. Integrating over these distribu- 
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tions and normalizing to N = N(< a-i) oc a 2 yields the total 
rate at which encounters are detected with an ELT 



r(a = 3/2) = 



dN(e,a) 



o j ai deda 
\5it^/2vq ( Gm; 



2o\ 



Tj(e,a)de da 

GM hh N 



5v„ 



a\ +«2+ \/ fl i fl 2 



(28) 



where a\ and a-i define the range in semi-major axis that is 
accessible to observations. 

Given the above expression for the encounter rate for a = 
3/2, we rely on scaling relations to estimate the rate for dif- 
ferent a. Since the encounter rate is proportional to the stellar 
density, T(a) ~ r(3/2)(r,,/r) Q - 3 / 2 . Thus, if a = 7/4, the rate 
of encounters is ~ 3 times larger than for a = 3/2. The time 
scale for detectable encounters is therefore 



'0.3 yr 



niiVo 



N 

Too 



2 x 10 5 M Q pc" 

V 02 ) 

V3000AU/ 



\0Mr. 



Oh 



100 km s" 



0.5 mas / \ 10 yr 



100 



-3/2 

(29) 

where we use the results of §0 that /V w 100, a\ ~ 200 AU, 
and «2 — 3000 AU, and have also assumed that the mass den- 
sity of background particles m, vo is constant and independent 
of mi. 

Therefore, assuming a density cusp dominated by ~ 1OM 
black holes, ~ 30 nearby stellar encounters will be detectable 
during ten years of monitoring with an ELT with 66 = 0.5 mas. 
Measurement of the frequency of detectable orbital deflec- 
tions r oc m.j is a direct test of the average mass of the dark 
remnants that proba bly dominate the mass density near the 
black hole jMorrisll993l [Miralda-Escud e & Gouldl2000l) but 
are otherwise not directly detectable. Since /V(< a) oc a 3 / 2 , 
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then by equation d28l >. T oc a 1 / 2 , i.e., the encounter rate in- 
creases with distance from the massive black hole. The stars 
at a > 3000 AU with detectable linear proper motion may 
therefore yield the strongest constraint on the mass function 
of stellar remnants, despite being below the threshold for de- 
tecting accelerated motion due to the massive black hole. 

4. METHOD 

In this section we describe how we generate mock ELT 
orbital data. We also describe our implementation of the 
Markov Chain Monte Carlo (MCMC) method which we use 
to estimate the uncertainties in the orbital parameters and map 
the shape of the likelihood surface. 

4.1. Parameter Estimation 

We are interested in estimating the uncertainties in the pa- 
rameters given proper motion and radial velocity information 
for a sample of N stars orbiting the massive black hole at the 
GC. Each star's projected orbit is described by six phase space 
parameters. The black hole mass, its 3-dimensional position, 
and the normalization and slope of the extended matter dis- 
tribution, contribute an additional six parameters. The en- 
tire parameter space of our model therefore has dimension 
J = 6N+6. 

Parameter estimation on a /-dimensional grid is not practi- 
cal. Since the computational cost of the grid-based approach 
increases exponentially with J, the parameter space becomes 
prohibitively large for even just two or three stars. By con- 
trast, the cost of the MCMC method scales almost linearly 
with J. 

We now briefly describe the basic ideas of the MCMC 
method and our choice of implementation. A general discus- 
sion of th e theory and ap plication of the MCMC a pproach 
is given in Gilks, Richardson, & Spiegelhalter ( 1996). Read- 
ers not interested in the d etails of our parameter estimation 
scheme can skip ahead to § 14.21 

Let D denote the observed data, 9 the model parameters, 
P(9) the prior distribution (which is uniform here), and L(D \ 9) 
the likelihood of detecting the data for a given set of parameter 
values. By Bayes's theorem the distribution of 9 conditioned 
on D is given by 



ir(9\D) = 



P(0)L(D\O) 
J P(9)L(D\9)d9 7 



(30) 



and is called the posterior distribution of 9. The statistical 
properties of the parameters such as means, moments and con- 
fidence contour levels, are entirely specified by ir{9\D). 

Explicit evaluation of the integral in the denominator of 
equation d30i is not practical in large dimensional models. 
The MCMC method avoids evaluating the integral by instead 
generating a Markov chain of parameter points the distribu- 
tion of which converges to the posterior distribution 7r(6>|D). 
The Markov aspect refers to the property that the probability 
distribution of the nth state (i.e., point) in the chain 9„ de- 
pends only on the previous state 9„ _i . It can be shown (e.g., 
iGilks. Richardson. & Spiegelhalteiti 19961) that the density of 
points in a Markov chain converges to ir(9\D) if the follow- 
ing criteria are satisfied: (1) the chain is irreducible, namely 
from any starting state 6o the chain can reach any non-empty 
set with positive probability in some finite number of itera- 
tions; (2) the chain is aperiodic in that it does not oscillate 
between different sets of states in a regular periodic fashion; 
(3) the chain is positive recurrent, meaning that if the initial 



value 6*o is sampled from the posterior then the expected time 
(i.e., number of iterations) to return arbitrarily close to state 6*o 
is finite. There are several algorithms for generating Markov 
chains that sa tisfy the above pr operties. We use the Metropo- 
lis algorithm (Metro polis et all!953l) in our numerical calcu- 
lations. 

Our implementation of the Metropolis algorithm is as fol- 
lows. 

1. Start a chain at t = with some initial state 9q. 

2. Generate a trial state 9' according to the jump proposal 
distribution q{9'\9,) (see below). Compute 



a{9 tl 9') = min 



UD\9') 
L(D\9 t ) 



(31) 



3. Sample a uniform random variable U that lies between 
(0,1). 

4. If U < a(9t,9') then set 9, +i =9' (i.e., accept the jump). 
If U > a{9,,9') then set 9, +{ = 9,. 

5. Increment t. 

6. Go to step #2. 

If the observational errors follow a normal distribution, 
L(D\9) ~ exp[-x 2 (6>,D)/2]. The x 2 (9, D ) statistic for a sin- 
gle star is given by 

2m m ^ [ [ Xi (9)- Xi (D)] 2 [ yi (9)-yi(D)] 2 \ 
X (0,D) = ^i -s- + 



K lvj(d)-Vj(D)] 2 



-E 



(32) 



where (x,y) is the astrometric position of the star, v its radial 
velocity, and a the corresponding measurement errors (i.e., 
&x,y = 30 ar, d f v = $ v )- We simultaneously fit to multiple stars 
by summing each star's % 2 to form a cumulative x 2 f° r the 
model. 

The jump proposal distribution q(9'\9 t ) is the probability 
of selecting a trial state 9' given the current state 9,. For the 
Metropolis algorithm one considers only symmetric proposals 
of the form q(9'\9,) = q(9,\9'). We choose to model the jump 
distribution as a multivariate normal distribution with mean 6, 
and constant covariance matrix C. 

Although the distribution of points in a chain is indepen- 
dent of the form of the jump distribution once the Markov 
chain has converged, the time it takes a chain to converge is 
sensitive to the jump distribution. To ensure an efficient run 
one must carefully chose the shape and step size of the jump 
distribution. An ideal jump distribution has a shape and step 
size that not only minimizes the convergence time but also 
samples the entire posterior distribution efficiently. In our im- 
plementation the shape of the jump distribution is determined 
by C and the step size is determined by a constant scale factor 
multiplying C. 

C is chosen such that the shape of the jump distribution is 
similar to that of the posterior distribution, although we again 
emphasize that the shape is only important for the efficiency 
of convergence. This ensures that the chain mixes well even 
in regions of degeneracy. To this aim, we compute the co- 
variance matrix that describes the shape of the x 2 surface in 
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the neighborhood of its minimum. We first compute the ap- 
proximate best fit parameter state 6*bf by minimizing \ 2 - We 
then specify a pilot covariance matrix C p that is purely di- 
agonal with variances given by a reasonable guess of the 1<t 
uncertainties for individual parameters. We draw a number 
(~ 1000) of pilot points from a multivariate normal distribu- 
tion with mean f?bf and covariance C p . Since the pilot points 
are within ~ 1 a of the \ 2 minimum, the shape of the \ 2 sur- 
face in the region of the points is approximately quadratic. We 
solve the linear least-squares problem by fitting a quadratic \ 2 
model to the points and obtain the approximate Fisher matrix 
that describes the curvature of the x 2 surface. We then deter- 
mine the eigenvalues and eigenvectors of the Fisher matrix. If 
any of the eigenvalues are negative, indicating that the shape 
of the surface is unconstrained in some direction, we generate 
more pilot points and redo the linear least-squares fit. Finally, 
we invert the resulting Fisher matrix to obtain the covariance 
matrix C. 

The constant scale factor that determines the step sizes must 
also be carefully chosen to ensure that the chain effectively 
explores the parameter space. If the steps are too small the 
chain does not mix well as it stays in one region of the pa- 
rameter space for long periods of time. If the steps are too 
large, the trial states are rejected frequently. For a multivari- 
ate normal jump distribution the most efficient step sizes are 
thos e for which ~ 25% of the jump proposals are accepted 
(see Gelman Hl995l) . We chose the (constant) jump scale fac- 
tor to optimize the acceptance rate. 

The first steps in a chain may be sensitive to the start- 
ing state 9() and are therefore not sampled from the pos- 
terior distribution. We discard these initial "burn in" 
points. To ensure that a chain has converged and is 
sampling the full posterior distribution we run multi- 
ple chains each starting at widely dispersed states. We 
tested for convergence with the G elman-Rubin test statistic 
JGilks. Richardson. & Spiegelhalterll996l) . 

4.2. Mock Data 

To generate a realistic set of orbital data we must determine: 
(i) the number of stars N we can detect and monitor with an 
ELT, (ii) the spatial distribution of these stars, (iii) the number 
of observations per year per star and, (iv) the observational 
errors in the stellar positions and velocities. In §[5]we showed 
that with a 30 meter ELT the position of the stars can be cen- 
troided to an astrometric precision 59io between 0.1-0.5 mas 
and the radial velocities measured to accuracies Sv^o between 
1-10 km s _I . We found that with such an ELT we can detect 
the accelerated proper motion of approximately 100 stars. We 
estimate that an integral-field spectrograph on an ELT enables 
a dedicated GC observing program to obtain position and ve- 
locities of each of the 100 stars roughly ten times per year. 

A realistic mock data set might therefore consist ofN = 100 
stars, observed over a ten year baseline with ten observations 
per year per star, with position and velocity measurements 
for each star accurate to 0.5 mas and lOkms -1 . Unfortu- 
nately, running our MCMC simulation on such a large data 
set was not feasible due to limits in computational speed. A 
typical run requires ~ 10 7 iterations (i.e., jumps) in order to 
fully sample the posterior distribution. This corresponds to a 
minimum of ~ 3 days on a desktop machine for just 20 stars 
(J ~ 126) with 100 points per star; a simulation with 100 stars 
takes approximately five times longer. However, one can ob- 
tain realistic results from a reduced sample size by properly 



TABLE 1 
Stellar Orbital Parameters 
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NOTE. — The listed parameters are: or- 
bital period (P), semi-major axis (a), eccentricity 
(e), pericenter distance (r m ; n ), and inclination (/). 
They are not all independent variables. 
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FIG. 2. — Astrometric positions of the 20 synthesized orbits to which we 
fit the model. The thick-lined portion of each orbit is the proper motion over 
the fitted 10 year baseline assuming 10 observations per year. 



scaling the \ 2 values (see equation I32> to emulate the full 
sample size. In particular, we construct a mock data set with 
N = 20 and multiply the \ 2 of each star by a factor of five. 

This approach yields realistic estimates of parameter uncer- 
tainties as long as the mock data set with N = 20 stars fairly 
represents the full data set with 100 stars. We minimize the 
effects of sample variance as follows. We first generate data 
for 1000 synthetic orbits. These orbits are drawn from the 
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distribution function of the power-law density profile assum- 
ing randomly oriented orbits and considering only those orbits 
with semi-major axes in the range detectable with an ELT (see 
§0. We generate mock data for these orbits assuming Gaus- 
sian position and velocity errors with dispersions 56 and Sv 
and a specific input model for the potential (e.g., black hole 
plus extended matter). For each individual star, we compute 
the difference in x 2 between the best-fit model (essentially 
the model used to generate the data) and the null hypothe- 
sis model (e.g., no extended matter). We then rank the stars 
by the size of this x 2 difference. We bin the 1000 stars into N 
bins according to their rank and randomly select one star from 
each bin. The resulting N stars form the set of orbits to which 
we fit. 

Table [2 lists the orbital parameters for one realization of 
a sample of 20 stars to which we fit. Given the orbital pa- 
rameters, we generate mock data by solving the equation 
of motion for each star (see equation In Figure [2] we 
show the astrometric positions of the 20 stars over the ten 
year observational baseline with ten epochs per year. The 
values of the input model parameters describing the poten- 
tial are: M bh = 4 x 1O 6 M , R = $ kpc, (x bh ,y bh ) = (0,0), 
M ext (r < 0.01 pc) = 6000M Q and either 7 = 1 .5 or 2. 

To test that the parameter uncertainty estimates are not af- 
fected by sample variance we ran simulations on several dif- 
ferent draws of 20 stars. As we show in §[5] the parameter 
uncertainties obtained are similar amongst the different data 
sets, suggesting that sample variance does not affect the re- 
sults. Thus, given the current uncertainties in an ELT's ulti- 
mate capabilities as well as the uncertainty in the exact nature 
of the stellar distribution at the GC, we conclude that to a rea- 
sonable approximation a mock data set comprised of N = 20 
stars with \ 2 values increased five-fold yields parameter un- 
certainties similar to that expected with observations by a 30 
meter ELT. 

We also show in the next section that the orbital parameter 
constraints scale with the measurement errors a and number 
of stars N as a/N 1 ^ 2 . Results for a wide range of assumed 
ELT capabilities (i.e., different 56 and Sv, different aperture, 
etc.) can therefore be computed by scaling the results of our 
fiducial 30 meter ELT model, using the relations between N, 
aperture D, 56, and Sv, given in § 

5. RESULTS 

In this section we investigate how well observations with 
an ELT constrain the structure of the GC. Our model of the 
GC and the orbits was described in § [3] We draw stellar or- 
bital parameters from a phase-space distribution determined 
by the model and use these orbits to synthesize mock ELT 
data (see §0). We then fit a model to the mock data and cal- 
culate the uncertainties i n th e parameters using the MCMC 
technique discussed in § 14.11 We show results for a 30 me- 
ter ELT with (56, Sv) = (0.5 mas, 10 km s" 1 ) and {56, 5v) = 
(0.1 mas, 2 km s -1 ). However, since the parameter uncertain- 
ties scale with measurement error a and number of monitored 
stars N as a/N 1 ' 2 , the results can be used to describe the capa- 
bilities of an ELT with different specifications. For exa mple , 
a 100 meter ELT will detect ~ 10 x as many stars (§ I2.2i : 
if the telescope has astrometric and spectroscopic errors that 
are smaller than those of a 30 meter telescope by a factor of 
five the parameter uncertainties will be ~ 10 x smaller. In this 
section, we estimate the limits that can be placed on the pa- 
rameters associated with the black hole including M b h and Ro 
(§ 15. 1> . as well as on the extended distribution of (dark) mat- 
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FIG. 3. — The constraints on Mbh and Rq obtainable with a an ELT as- 
suming an extended matter distribution with 7 = 1.5 (results are similar for 
7 = 2). The solid contours show the 68%, 95%, and 99.7% confidence lev- 
els assuming an astrometric limit of 59 = 0.5 mas and a spectroscopic limit 
of Sv = 10 km s for the draw of 20 stars shown in Table fTI The line-dot 
contour shows the 99.7% confidence level for a different draw of 20 stars. 
The dashed contour shows the 99.7% confidence level for smaller astromet- 
ric and spectroscopic limits of 59 = 0. 1 mas and 5v — 2 km s - ' . An ELT will 
constrain both M(,h and Rq to better than 0.1%. 
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FIG. 4. — The constraint on the extended matter distribution obtainable 
with an ELT. Shown are the 68%, 95%, and 99.7% confidence levels on the 
enclosed mass and slope of an extended matter distribution assuming an as- 
trometric limit of 59 = 0.5 mas and a spectroscopic limit of 5v = 10 km s . 
The input models have power-law slope of 7 = 1.5 and 7 = 2 and an input en- 
closed mass of 6000Mq within 0.01 pc. The dashed contour is the constraint 
at the 99.7% level for measurement errors that are a factor of five smaller. 



ter near the black hole (§ 15. 21 . We discuss the dependence 
of the limits on the astrometric and spectroscopic precision 
of the observations. We also investigate whether relativistic 
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FIG. 5. — An ELT's sensitivity to post-Newtonian effects assuming an as- 
trometric limit of 89 = 0.5 mas and a spectroscopic limit of 8v = 10 km s . 
Shown is the uncertainty in the speed of light and the extended matter mass 
as obtained by including post-Newtonian corrections to the equations of mo- 
tion. The Roemer delay and special relativistic effects are not included in the 
model in order to demonstrate that general relativistic effects of order (v/c) , 
including the prograde precession, are detectable with an ELT. The line styles 
are the same as in Figurel3l 

corrections to the Keplerian motion can be detected at the GC 

(§122). 

5.1. Measuring Mbh and Rq 

In Figure [5] we show the constraints an ELT will place on 
Mbh and Rq. For an astrometric limit of 5 6 = .5 mas and a 
spectroscopic limit of 5v = 10 km s~ (see § 12. U the fractional 
uncertainties in Mbh and Rq are less than 0.1% at the 99.7% 
level. This is a factor of ~ 100 times better than present uncer- 
tainties. The result is robust in that simulations with distinct 
mock data sets of 20 stars, drawn in the fashion described in 
§ 14.21 produce similar uncertainties in the model parameters. 

For astrometric and spectroscopic limits that are a factor 
of five smaller the fractional uncertainties in Mbh an d Ro are 
smaller by almost a factor of five. The uncertainties in Mbh 
and Rq scale almost linearly with the measurement errors for 
observations at this precision. We also verified that the uncer- 
tainties scale with N as roughly N~ x l 2 . 

Observations with a 30 meter ELT will therefore constrain 
the distance to the GC to within a few parsecs and the mass of 
the black hole to within a few thousand solar masses. We dis- 
cuss the implications of measuring Rq to such high accuracy 
in SB 

5.2. Measuring the Extended Matter Distribution 

In Figure@]we show the constraints an ELT will place on the 
extended matter distribution for input power-law models with 
M ext (r < 0.01 pc) = 6000M Q and 7 = 1 .5 or 7 = 2. We chose 
these distributions in order to conform to the extrapolation of 
the observed stellar density distribution and to theoretical es- 
timates of dark matter clustering (see §[Q. We find that for an 



ELT with 56 = 0.5 mas and 5v=\0 km s~ one can detect such 
extended matter distributions, yielding measurements of M ext 
and 7 that are accurate to 20-30% (i.e., <5M ext ~ 1500M Q 
and Sj ~ 0.5). Since the amplitude of the N ewto nian retro- 
grade precession varies linearly with M ext (§ 13. 21 . the frac- 
tional uncertainty is 5M ext /M ext oc £A</> Newt / A0 Ne „ t oc M~ xl , 
where 5 A^Newt is set by the astrometric precision. Thus 5M ext 
is independent of M ext so that an extended matter distribution 
is detectable (i.e., observations yield a lower bound) for 56 = 
0.5 mas and 5v = 10 km s -1 as long as M ext (r < 0.01 pc) > 
5M ext ~ 1500M Q . Such an ELT will therefore place interest- 
ing constraints on the extended matter at the GC. 

5.3. Measuring Relativistic Effects 

As discussed in § 13.31 order of magnitude estimates suggest 
that post-Newtonian corrections to the equations of motion, 
involving terms of order (v/c) 2 , are measurable with an ELT 
with astrometric resolution of 59 < 0.5 mas. In an effort to 
demonstrate this more quantitatively, we allow the speed of 
light to be a parameter in our model and examine how well we 
recover its value. We purposely do not include relativistic cor- 
rections to the observed motion associated with propagation 
effects (e.g., the Roemer time delay and other higher-order 
corrections) so that we can examine the detectability of (v/c) 2 
general relativistic corrections to the orbital dynamics such as 
the prograde precession of the major axis position. In Figure 
13 we show the constraint on c as a function of M enc . Post- 
Newtonian effects are observable, as c is measured to ~ 5% 
accuracy. Since v/c < 0.2 for all stars in the sample fS I3.lt 
the few percent constraint on c suggests that while the (v/c) 2 
effects are measurable, the (v/c) 3 effects are not. The orbital 
precession due to black hole spin is of order (v/c) 3 (§ 13.41 
and detecting it with an ELT with 59 = 0.5 mas requires the 
favorable discovery of a star on a compact and highly eccen- 
tric orbit. Based on estimates of the signal-to-noise from a 
spin-induced apocenter shift (equation 1151 ). an astrometric 
precision of ~ 0.05 mas is needed to reliably detect the black 
hole spin. 

The degeneracy between c and M enc is a consequence of 
the degeneracy between the prograde relativistic precession 
and the retrograde Newtonian precession. Decreasing c in- 
creases the amount of prograde motion A</> pro while increas- 
ing M enc increases the amount of retrograde motion A^Newt- 
The two effects compensate for one another over a range of 
c and M enc . The degeneracy is broken at sufficiently extreme 
values of M enc because the relativistic and Newtonian effects 
each induce a distinct precessional shape. 

6. CONSTRAINTS ON GALACTIC STRUCTURE FROM 
MEASUREMENTS OF R a 

The distance to the GC, Rq, is a funda mental parameter in 
model s of the Milky Way structure. As Oiling & Merrifieldl 
(2001) note, models of the Milky Way exhibit strong inter- 
relations between the Galactic constants (Rq and the local 
Galactic rotation speed Oo), the shortest-to-longest axis ra- 
tio, q = c/a, of the dark matter halo, and the local stellar col- 
umn density £*. The determination of q is of particular inter- 
est since different models of dark-matter and structure forma- 
tion scenarios predict different values for q. Cold dark-matter 
simulations t ypically produce galactic halos that are triaxial 
JWarren etalJl!992t Uin g & Sutol 12002) although these tend 
become oblate under the influence of the d issipative infall 
of gas resulting in halos with q ~ 0.5 (Dubinski 1994J). Al- 
ternatively, hot dark-matter models predict round halos with 
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q ~ 0.8 (Peebles 1993) wh ile some baryonic dark m atter mod- 
els imply q ~ 0.2 (Pfenniger. Combes. & Martinet 1994). As 
we now discuss, determining Rq to 0.1% via monitoring of 
stellar orbits at the GC with an ELT enables an extremely pre- 
ci se measurement of q in the Milky Way. 

Oiling & Merrifield (2000) demonstrate that there is signif- 
icant uncertainty in existing estimates of q in galaxies due to 
both the limited amount of data available for measuring q and 
the fact that different measurement techniques have yielded 
systematically different values. Presently, the situation is not 
any better for our own Galaxy, with plausible values lying in 
the range 0.3 < q < 1. 

The measurement of q in the Milky Way entails measur- 
ing the Galaxy's radial mass distributio n and the degree to 
which this mass distribution is flattened. Oiling & Merrifield! 
(2000) show that the uncertainty in q in the Milky Way is al- 
most entirely due to the l arge errors in the Gal actic constants 
6 and R . Indeed, lolling & Merrifield! (120011) show that the 
fractional uncertainty in q is nearly twice the fractional uncer- 
tainty in Oq. Therefore, a precision measurement of the Sun's 
proper motion with respect to the GC in combination with a 
precision measu rement of Rq tightly constrain s 8o and hence 
q. According to lSalim. Gould. & Olling l J2002) future astro- 
metric surveys will be able to measure the Sun's proper mo- 
tion \x = V /Rq to within several microarcseconds, correspond- 
ing to 0.1% accuracy. Here V = 9o + V© ~ 220 km s -1 is the 
sum of the rotation speed of the local standard of rest and the 
Sun's motion relative to it. The uncertainty in Oo will be the 
dominant error in V; V© is already known to an accuracy of 
0.6 km s -1 from the Hipparcos catalogue ( Dehnen & Binnev 
(19981) . Thus, the monitoring of stellar orbits at the GC with 
an ELT in conjunction with future astrometric survey missions 
will constrain the Milky Way's dark matter halo shape param- 
eter q to a few tenths of a percent. 

7. CONCLUSIONS 

We have examined a variety of experiments that can be 
achieved through the infrared monitoring with an ELT of stars 
within a few thousand AU of the GC. The astrometric limit of 
a 30 meter ELT is conservatively 0.5 mas and possibly as high 
as 0.1 mas. By comparison, the astrometric limit of current 
observations is 1 - 2 mas. 

The greater point-source sensitivity and spectral resolution 
of an ELT enables the measurement of radial velocities with 
errors < 10 km s -1 . At present, of the ~ 10 stars with mea- 
sured accelerated proper motions, spectral lines have been 
detected only in SO-2, with radial velocity uncertainties of 
~ 30 km s -1 . Measuring the radial velocities of stars breaks 
the degeneracy between mass and distance and thus yields a 
direct measurement of the distance to the GC. If the spectra 
of fainter stars can be obtained, the detection of deep molecu- 
lar lines will improve upon the velocity estimates by an addi- 
tional factor x 10. The solar type stars that will be detectable 
with an ELT may therefore yield radial velocity uncertainties 
considerably smaller than 10 km s -1 . 

A 30 meter ELT will be able to detect stars down to a K- 
band magnitude of K ~ 22, approximately four magnitudes 
fainter than currently possible. Due to confusion, it will be 
difficult to detect still fainter stars. Using measurements of 
the A"-band luminosity function within the inner 1" of the 
GC, we estimate that such an ELT will detect the acceler- 
ated motion of ~ 100 stars with semi-major axes in the range 
200 < a < 3000 AU. Current observations are limited to the 



detection of — 10 stars, all with a > 1000 AU. We find that 
the number of stars with detectable accelerated motion scales 
with the aperture of an ELT as N ~ 100(£>/30 m) 2 . 

Given the observational capabilities of an ELT and the 
likely, albeit at low masses largely uncertain, stellar environ- 
ment at the GC, we constructed a plausible sample of stellar 
orbits. The model includes the dynamical contribution of an 
extended distribution of dark matter around the black hole that 
is composed of stellar remnants and CDM. We find that for 
measurements at the precision obtainable with an ELT the un- 
certainty in the model parameters scale with the measurement 
errors a (i.e., 56, 5v) and the number of monitored stars N 
as roug hly cr/JV 1 / 2 . Thus, while we focus on the capabilities 
of a diffraction limited 30 meter ELT with 56 = 0.5 mas and 
5v = 10 km s _I , our results can be used to determine the capa- 
bilities of an ELT with different specifications. For example, 
a 100 meter ELT will detect ~ 10 x as many stars so that if 
it has astrometric and spectroscopic errors that are smaller by 
a factor of five, the measurement accuracy in the parameters 
will improve by a factor of approximately ten. 

We find that with a 30 meter ELT the parameters Mbh and 
Ro will be measured to an accuracy better than 0.1%. Deter- 
mining Ro to within a few parsecs will significantly constrain 
models of the Galactic structure as it aids the precise measure- 
ment of the dark matter halo shape. 

While current observations of stellar proper motions are 
compatible with Keplerian motion, a number of dynamical 
effects produce significant deviations, including the Newto- 
nian retrograde precession, the relativistic prograde preces- 
sion, frame dragging due to the black hole spin, and inter- 
stellar interactions involving nearby encounters. All but the 
frame dragging effect produce non-Keplerian motions that 
are detectable with a 30 meter ELT. Unfortunately, the spin 
of the massive black hole at the GC will probably be out of 
reach to kinematic studies unless an astrometric precision of 
~ 0.05 mas is achieved. 

The presence of an extended distribution of matter results 
in a Newtonian retrograde precession due to differences in 
the amount of mass enclosed within an orbit's pericenter and 
apocenter. We considered extended matter density profiles 
consistent with current observations of the stellar distribution 
at the GC. We modeled the distribution as a power-law profile 
normalized such that M en (r < 0.01 pc) = 6000M© and with 
slope 7 = 1 .5 or 2. Standard models of dark matter clustering 
about a massive black hole predict similar profiles. An orbit 
monitoring program with a 30 meter ELT will constrain the 
mass and slope of such profiles to ~ 30% accuracy. Thus, 
monitoring orbits with an ELT provides a probe of the ex- 
tended matter distribution within ~ 10 4 Schwarzschild radii 
of the massive black hole at the GC. 

We also calculated the rate at which the monitored stars 
experience detectable deflections due to stellar gravitational 
scattering encounters with background compact remnants. 
We considered a detection threshold set by the minimum de- 
tectable change in the velocity of a monitored star. For a den- 
sity cusp dominated by ~ 10M Q black holes, ~ 30 nearby 
stellar encounters will be detected by a 30 meter ELT over a 
ten year observing baseline. This will confirm the presence 
of a cusp of compact remnants at the GC and enable the mea- 
surement of the remnants' masses. 
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